BOUNDARY CLUSTERED LAYERS NEAR THE HIGHER 
CRITICAL EXPONENTS 



NILS ACKERMANN, MONICA CLAPP, AND ANGELA PISTOIA 

Abstract. We consider the supercritical problem 

—Am = \u\ p ~ 2 u in Q, u = on dQ, 
where Q is a bounded smooth domain in M. N and p smaller than the crit- 
ical exponent 2* N k := %^zirzl for the Sobolev embedding of H 1 (M N ~ k ) in 
L q (S. N ~ k ), 1 < k < N — 3. We show that in some suitable domains Q there 
are positive and sign changing solutions with positive and negative layers which 
concentrate along one or several fe-dimensional submanifolds of dQ as p ap- 
proaches 2* N , from below. 
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1. Introduction 
Consider the classical Lane-Emden-Fowler problem 
(1) Av + \v\ p - 2 v = in V, v = ondV, 

where T> is a bounded smooth domain in WL N and p > 2. 

It is well known that when p is smaller than the critical Sobolev exponent 2* := 
-^32 ; compactness of the Sobolev embedding ensures the existence of at least one 
positive solution and infinitely many sign changing solutions. In contrast, existence 
of solutions to problem ([1]) when p > 2* is a delicate issue. Pohozhaev's identity [22] 
implies that problem ([T]) does not have a nontrivial solution if the domain T> is 
strictly starshaped. On the other hand, Kazdan and Warner showed in [T3] that if 
the domain T> is an annulus, problem (JTJ) has infinitely many radial solutions. 

For the critical case p = 2* Bahri and Coron [T] proved that a positive solution 
of ([1]) exists if the domain T> has nontrivial reduced homology with Z/2-coefficients. 
Moreover, it was proved by Ge, Musso and Pistoia [TT] and Musso and Pistoia [TB] 
that, if V has a small hole, problem (JlJ has many sign changing solutions, whose 
number increases as the diameter of the hole decreases. Multiplicity results are 
also available for domains which are not small perturbations of a given domain, but 
have enough, possibly finite, symmetries, as proved by Clapp and Pacella [8] and 
Clapp and Faya [6]. 

The almost critical case p — 2* ± e, with e positive and small enough, has 
been widely studied. The slightly subcritical case p = 2* — e was considered by 
Bahri, Li and Rey [2] and Rey [23], who showed the existence of positive solutions 
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which blow-up at one or more points of T> as e — > 0. A large number of sign 
changing solutions with simple or multiple positive and negative blow-up points 
were constructed by Bartsch, Micheletti and Pistoia 3 , Musso and Pistoia |17j . 
and Pistoia and Weth [21]. For the slightly supercritical case p = 2* + e existence 
and nonexistence of positive solutions with one or more blow-up points has been 
established by Ben Ayed, El Mehdi, Grossi and Rey [9], Pistoia and Rey [20], and 
del Pino, Felmer and Musso [5]. 

Unlike the critical case, in the supercritical case p > 2* the existence of a non- 
trivial homology class in T> does not guarantee the existence of a nontrivial solution 
to |T]). In fact, for each integer k such that 1 < k < N — 3, Passaseo [181119] 
exhibited a bounded smooth domain in R , homotopically equivalent to the k- 
dimcnsional sphere, in which problem (p} does not have a nontrivial solution for 
P > 2*n k : ~ iv-fc~-2 ■ Note that 2^, fe is the critical Sobolev exponent in dimension 
N — k. Examples of domains with richer homology were recently given by Clapp, 
Faya and Pistoia [7] , where it was shown that for p > 2* N k there are bounded 
smooth domains in M. N whose cup- length is k + 1, in which problem (fTJ) does not 
have a nontrivial solution. On the other hand, for p = 2* N k existence of infinitely 
many solutions in some domains has been recently established by Wei and Yan [25] . 
Further multiplicity results may be found in [7]. 

In [TU] del Pino, Musso and Pacard considered the case p = 2* N 1 — e and proved 
that for some suitable domains V, if e is positive, small enough and different from 
an explicit set of values, problem (fT} has a positive solution which concentrates 
along a 1-dimensional submanifold of the boundary dT>. In the same paper, the 
authors ask the question whether one can find solutions which concentrate at a 
/s-dimensional submanifold for p slightly below 2^, k . More precisely, they ask the 
following: 

Problem 1.1. Given 1 < k < N — 3, are there domains T> in which problem (QJ) 
has a positive solution v p for each p < 2^- k with the property that these solutions 
concentrate along a fc-dimensional submanifold of the boundary dT> as p — > 2* N k ? 

Having in mind that when p approaches the first critical exponent 2* from below 
a large number of sign changing solutions exist, another question arises naturally: 

Problem 1.2. Given 1 < k < N — 3, are there domains T> in which problem (|7J) has 
a sign changing solution v p for each p < 2^ k with the property that these solutions 
concentrate along a fc-dimensional submanifold of the boundary dT> as p — > 2* N > ? 

In this paper, we give a positive answer to both questions. In particular, for 
each set of positive integers k%, . . . , k m with k := k\ H — • + k m < N — 3 we exhibit 
domains T> in which problem (fTJ) has a positive solution for each p < 2^ k and, 
as p 2* N k , these solutions concentrate along a fc-dimensional submanifold M of 
the boundary &D which is diffcomorphic to the product of spheres S fcl x • • • x E> km . 
Moreover, problem (fT]) has also a sign changing solution with a positive and a 
negative layer, both of which concentrate along M as p — > 2* N k . This follows from 
our main results, which we next state. 

Fix k\, . . ., k m € N with k := k\ + • • • + k m < N — 3 and a bounded smooth 
domain fl in M. N ~ k such that 

(2) He {(xi,...,x m ,x') € R m x R N - k - m : Xl >0, i = l,...,m}. 



BOUNDARY CLUSTERED LAYERS NEAR THE HIGHER CRITICAL EXPONENTS 3 



Set 
(3) 

V := {(tfV..,i/ TO ,«) g K fcl+1 x ••• x R k ™ +1 x R^-k- m . fly 1 ] >J8 ) g ft}. 

I? is a bounded smooth domain in which is invariant under the action of the 
group T := 0{ki + 1) x • • • x 0(k m + 1) on R N given by 

for every g f g 0(fc 4 + 1), y i g K fc ' +1 , z g l"-^™. Here, as usual, O(d) denotes the 
group of all linear isometries of M. d . For p = 2* N , — e we shall look for T- invariant 
solutions to problem ((TJ), i.e. solutions v of the form 

(4) v(y 1 ,...,y m ,z) = u(\y 1 \,...,\y m \,z). 

A simple calculation shows that v solves problem ([I]) if and only if u solves 

~Au-jr— = \u\ p - 2 u in SI, u = ondn. 

This problem can be rewritten as 

— div(a(x)Vu) = a(x)\u\ p ~ 2 u in SI, u = on dtt, 

where a{x\, . . . , xjv— fc) := ^i 1 •••2 ; m™- Note that 2^ fc is the critical exponent in 
dimension n := N — k which is the dimension of SI. 

Thus, we are lead to study the more general almost critical problem 

(5) — div(a(x)Vw) = a(x) \u\ n ~ 2 u in SI, u = on dfl, 

where SI is a bounded smooth domain in R™, n > 3, e is a positive parameter, and 
a g C 2 (S1) is strictly positive on SI. 

This is a subcritical problem, so standard variational methods yield one positive 
and infinitely many sign changing solutions to problem ([5]) for every e g (0, 335)) 
cf. Proposition 4.1 in [7]. Our goal is to construct solutions u e with positive and 
negative bubbles which accumulate at some points £1, . . . of 9S1 as e — > 0. They 
correspond, via to T-invariant solutions u e of problem (fTJ) with positive and 
negative layers which accumulate along the /c-dimensional submanifolds 

Mj ■.= {(y 1 ,...,y m ,z)eR k ^ 1 x---xR k -+ 1 xR N - k - m : ((y 1 1 , . . . , \y m \ , z) = 

of the boundary of T> as e — > 0. Note that each Mj is diffeomorphic to § fel x • • • x E> km 
where S d is the unit sphere in R d+1 . 

We will assume one of the following conditions. 

(al) There exist k nondegenerate critical points £1 , . . . , £ K g <9S1 of the restriction 
of a to <9S1 such that 

(Va(&), >0 Vi = l,...,/s, 

where is the inward pointing unit normal to 917 at 
(a2) There exists a critical point £0 S 90 of the restriction of a to <9S1 such 
that (Va(£o), K£o)) > 0, and vectors t%, . . . , r n _i g M n such that the set 
{^(£0)1 r i7 • ■ • i Tn— 1} is orthonormal and SI and a are invariant with respect 
to the reflection gi on each of the hyperplanes £ + { T i = 0}, i.e. 

Qi(x) g S7 and a(gj(a;)) = o(.t) \/x g Q, 
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i = 1, ...,n — 1 , where 

ft(£o + (», ^) v+ (x,ri) n H h (x,Ti)Ti H h (x,T n _i) r„_i) 

= £o + (», f) v + (x,n) t\ H (x,Ti) n H h (a;,r n _i) r„_i. 

and f := f(£o) is the inward pointing unit normal to d£l at £o- 
For each 6 > 0, £ € R™, we consider the standard bubble 

I7 4|£ (a;) := [n(n-2)]^ — — . 

(5 2 + \x-^) — 

We prove the following results. 

Theorem 1.3. Assume that (al) holds true. Then there exists eo > such that, 
for each Ai,...,A re € {0,1} and e £ (0, eo), problem |5} has a solution u e which 
satisfies 

K 

u t (x) = ^(-l) A *t^ i ., fii .(a:) + o(l) m D^ 2 (n), 
i=i 

with 

e~5i=2(5 i)e -> dj > and £ l>e -> & £ <9fi, 
/or eac/i z = 1, . . . , k, as e — > 0. 

Theorem 1.4. Assume that (a2) holds true. Then there exists e > smc/i i/iai, 
for each e € (0,eo), problem ([5]) has a sign changing solution u e which is invariant 
with respect to each reflection Qi, i — 1, ...,n — 1, and satisfies 

lie fa) = f^i,«,^i, e (a:) - f^a.e^.eC 35 ) + c 1 ! 1 ) * n £)1:2 (°)) 
with i 

e~~<5 ii£ -> d, > 0, £ i>e = £o + e*i,eK£o) a«d *i,e -> *« > 0, 
for each i = 1, 2, as e — > 0. 



Theorem 11.41 states the existence of a sign changing solution whose two blow-up 
points (one positive and one negative) collapse to the same point £o of the boundary 
of O under the symmetry assumption (a2). 

Some interesting questions arise: 

Problem 1.5. Is it possible to find sign changing solutions with k > 3 blow-up 
points with alternating sign which collapse to the point £o? 

Problem 1.6. Is it possible to find a sign changing solution with one positive 
and one negative blow-up point which collapse to the point £o in the more general 
case when £o is a nondegenerate critical point of a constrained to dQ such that 
(Va(£o), ^(£o)) > 0> without any symmetry assumption? 

The reason for including the symmetry assumption (a2) in Theorem 11.41 is that 
it allows to simplify the computations considerably (see Remark 12.61) . 

In the following two theorems we assume we are given k\ , . . . , k m 6 N with 
k := &! + ••• + k m < N — 3 and a bounded smooth domain £1 in M. N ~ k which 
satisfies ([2]). We set a(xi, . . . ,XN-k) '■= xl 1 ■ ■ ■ x k ™ , V as in ([3]), p = 2^ fc — e, 
r := 0(h + 1) x • • • x 0(k m + 1) and 

U SA (y\...,y m ,z) ■.^U 5 ^\y 1 \,...,\y m \ 1 z) 

for 5 > 0, f G E. N ~ k . 
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Theorem 1.7. Assume that (al) holds true for a and fl as above. Then there 
exists €q > such that, for each Ai, . . . , X K G {0, 1} and e G (0, eo), problem ([T]) has 
a I '-invariant solution v e which satisfies 



with 



i=l 



e -- 2 (Si, e -> di > and &, e -> & G <9f2, 
/or eac/i i = 1, . . . , re, as e — > 0. 

Theorem 1.8. Assume that (a2) /io/ds irite /or a and 17 as above. Then there 
exists e > suc/i £/ia£, /or eac/i e G (Q,eo), problem |l| /ias a T-invariant sign 
changing solution v e which satisfies 

v e (x) = U Sl ,.,ti,. {x) - (a;) + o(l) m D 1,2 (T>), 

with 

— n_1 

e "- 2 S i:€ -> d,; > 0, £ i>e = £o + e*i,eK£o) and t it6 -> ^ > 0, 
/or eac/i i = 1,2, as e — > 0. 

By the previous discussion Theorems 11.71 and 11.81 follow immediately from The- 
orems 11.31 and 11.41 The proof of Theorems 11.31 and 11.41 relies on a very well known 
Ljapunov-Schmidt reduction procedure. We shall omit many details on this proce- 
dure because they can be found, up to some minor modifications, in the literature. 
We only compute what cannot be deduced from known results. 

The outline of the paper is as follows: In Section [5] we write the approximate 
solution, sketch the Ljapunov-Schmidt procedure and use it to prove Theorems 11.31 
and 11.41 In Appendix B we compute the rate of the error term and in Appendix C 
we estimate the reduced energy. In Appendix A we give some important estimates 
on the Green function close to the boundary. 

2. The variational setting 

We take 



(u, v) := / a(x)Vw • Vu dx, \\u\\ := ( / a(x) | Vit| 2 dx I 
Jn \Jn J 



1/2 



as the inner product in Hj(f2) and its corresponding norm. Since a is strictly 
positive and bounded in f2 they are well defined and equivalent to the standard 
ones. Similarly, for each r G [l,oo), 

l/r 

a(x) \u\ r dx 



is a norm in L r (fJ) which is equivalent to the standard one. 

Next, we rewrite problem (|S|) in a different way. Let i* : L^(fi) -> Hj(fi) be 
the adjoint operator to the embedding i : Hj(i7) <^-> L n ~ 2 (Q), i.e. i*(u) = v if and 
only if 

(v,ip) = / a(x)u(x)ip(x)dx for all ip G C^°(Q) 



if and only if 

— div(a(a;)Vu) = a(x)u in Q, v — on 8Q. 
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Clearly, there exists a positive constant c such that 

\\i*(u)\\ < C\\U\\ 2^ VttGL^F2(^). 

Setting p :— an d /e( s ) : = l s | P 2 e s , problem §5§ turns out to be equivalent to 
(6) « = i*(/e(«)), ueHj(fi). 

n-2 

Set /(s) := /o(s) and a n := [n(n — 2)] 4 . Let 



U s ,i~a n S>0, (el", 

be the positive solutions to the limit problem 

-Au = f(u), uei^QfT). 

Set 

:= ~W = a 2 { p + \x-W 

and, for each j = 1, . . . , n, 

^ W :A^„ ( „- 2 )^ (<3+ ^ : ^ r/2 . 

Recall that the space spanned by the (n + 1) functions is the set of solutions 
to the linearized problem 

-Aip = (p- l)U%~ 2 ip in R n . 

Let PW denote the projection of the function W £ D 1 ' 2 (W l ) onto Ho(fi), i.e. 

APW = AW in fl, PW = on dfl. 

We look for two different types of solutions to problem ([5]) . The solutions found 
in Theorem 11.31 are of the form 



(7) u t = Y J (-^PUs,^+^, 

i=l 

for fixed € {0, 1}, where the concentration parameters satisfy 

(8) = e^^di for some di > 0, 
and the concentration points satisfy 

(9) £i, e = Si + rjiv{si) where s; G dQ and 7}% = eti for some U > 0. 

Here and in the following v(si) denotes the inward unit normal to the boundary 
d£l at the point s$. 

On the other hand, the solutions found in Theorem 11.41 are of the form 

e 

(10) u e = ]T(-iy +1 p;y AWi , c + 0, 

j=i 

where the concentration parameters satisfy (|8|), while the concentration points are 
aligned on the line £ := {£o + r K£o) : r £ K}, namely 

(11) £i, € = £ + Vi^(^o) where m = eti for some < t\ < ■ ■ ■ < fy. 
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Next, we introduce the configuration space A where concentration parameters 
and concentration points lie. For solutions of type ([7]) we set s = (si, . . . , s K ) £ 
(dn) K , d = (di,. . .,d K ) G (0, +oo) K , and t = (ii,. . . ,t K ) £ (0, +oo) K , and so 

A := {(s, d, t) G (W) K x (0, +oo) K x (0, +oo) K : Si ^ Sj if i ^ j} , 

while for solutions of type (|10p. we fix s = (£q, • ■ • , £o) and we set d = (di, . . . , d^) G 
(0, +oo) £ , and t = (ti, . . . , tg) £ (0, +oo)*, and so 

A := {(d, t) G (0, +oo) e x (0, +oo) e : ti < • • • < t e } . 

In each of these cases we write 

K I 

V.AX ■■= ^(-l) A *PC/<5 i)£ . and V B , d ,t = Vd, t := ^(-l^PE/^ 

i=l i=l 

respectively. 

The rest term <fi belongs to a suitable space which we now define. For simplicity 
we write ty\ := ify ^ t. e with <5, ie as in ((8|) and £i i£ as in (j9|) or (fTTj) . 
For solutions of type © we introduce the spaces 

#s,d,t := span{P?/>- : i = 1, . . . , k, j = 0, 1, . . . , n}, 
#sUt == {<^> e Hj(O) : (0, P$) = 0, i = 1, . . . , k, j = 0, 1, . . . , n} . 

Note that for £j. e as in (jTTj) the functions P?/>f are invariant with respect to 
the reflections Qi given in (a2). So for solutions of type (|T0]) we define the space 
^s,d,t as above and -K^j t as the orthogonal complement of K s dt in the subspace 
of all functions in Hj(51) which are invariant with respect to Q\,...,Q n —\- Then 
we introduce the orthogonal projection operators n s d,t and IT^ d t in Hq(Q) with 
ranges K sd t and K^~ d t , respectively. 

As usual, our approach to solve problem © will be to find a (s, d, t) G A and a 
function <f> £ K^ d t such that 

(12) n^ t (v s , d<t + 4>-i* [f e (v Sld ,t + = o 

and 

(13) n s , d> t (v;,d,t + 4 - i* [/ e (K,d,t + 0)]) = o. 

First we shall find, for each (s, d, t) G A and small e, a function tj> £ K^~ d t such 
that (|I2I) holds. To this aim we define a linear operator £ S! d,t : i^j j t — * ^ s d t by 

Ls,d,t0:=^»-n^ d>t i* [/ / (7., d ,t)0] • 
The following statement holds true. 

Proposition 2.1. for any compact subset C of A £/iere exist eo > and c > 
smc/i t/iat /or each e £ (0, eo) and (s, d, t) G C i/ie operator i Si d,t *s invertible and 

||£s,d,t0|| >c||0|| V^e^ d ,. 

Proof. We argue as in Lemma 1.7 of [15] , □ 
Now we are in position to solve equation (|12[) . 
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Proposition 2.2. For any compact subset C of A there exist eo,c, c > such that 
for each e £ (0, eo) and (s, d, t) £ C there exists a unique (f>^ d t £ K^~ d t such that 
(|12[) holds and 

(14) ||# id)t || <«*+". 
Proof. We estimate the rate of the error term 

(15) #s,d,t := n s ^ d , t (V s , d , t - i* [/e(T4, d ,t)]) 

in Appendix B. Then we argue exactly as in Proposition 2.3 of [3]. □ 
The critical points of the energy functional J € : Hj(57) — > R defined by 



i / a(a;)|Vw| 2 ^ — / a(a:)|u| p - 6 da: 

n si 



are the solutions to problem ([5]) . We define the reduced energy functional J e : A 
R by 

J e (s,d,t) := J e {V s , d ,t + ^ jd ,t) 
The critical points of J e are the solutions to problem (fT3| . 



Proposition 2.3. TTie function V s .d,t + 0s d t * s a critical point of the functional 
J e if and only if the point (s, d, t) is a critical point of the function J e . 

Proof. We argue as in Proposition 1 of [2]. □ 

The problem is thus reduced to the search for critical points of J e , so it is 
necessary to compute the asymptotic expansion of J e . 

Proposition 2.4. In case it holds true that 

(16) 

K 

J e (s, d, t) = (ci + c 2 e log e) ^a(sj) 



i=l 



n-2 



+ o(e), 



^-uniformly on compact sets of A. Here the Ci's are constants and C4,C5,C6 are 
positive. 

Proof. The proof is postponed to Appendix C. □ 
Proposition 2.5. In case (jTQl) it holds true that 

(17) J e (s, d, t) = J e (d, t) = o(fo) [ci + c 2 e loge + c 3 e] + etf (d, t) + o(e), 
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C° -uniformly on compact sets of A. Here 



(18) *(d,t) :=C4<Va(£o),K£o))I> + C5a(£o): 



i=i 



E 



\u-t 3 \ n - 2 \ti + tj\ n -* 



c 6a(Co)^ lo S dl 



where the c$ 's are constants and 04,05, C6 are positive. 
Proof. The proof is postponed to Appendix C. 

Proof of Theorem \1.3[ Firstly, by Proposition 12.41 we get 



□ 



J £ (s, d,t) = (ci + c 2 e log e) ^a(sj) + 0(e), 



i=i 



C^-uniformly on compact sets of A. Then, since £1, • • • ,£« are non degenerate critical 
points of a constrained to the boundary of O, if e is small enough there exist 
s e := (si, e , . . . , s K)£ ) such that each s,- )f — ► as e goes to zero, and V s J e (s e , d, t) = 0. 
Secondly, by Proposition 12.41 we also get 



J £ (s e ,d,t) - (ci + c 2 e log e) ^a(s ije ) 



; E 



C3a(s»,e) + c 4 (Va(s ¥ ), i/(sj, e ))*i 



+ csa(si :£ ) ( — 



c 6 a(si !£ ) log 



+ o(e) 



= e E 



C3a(6) + c 4 (Va(&), v(£i)}ti 



c S a(6) ( 



n-2 



c 6 a(&)log 



•o(e). 



It is easy to verify that the function 



(d,t)->£; 



i=l 



c 4 (Va(&), f(^))ii + c 5 a(^) 



n-2 



c 6 a(£j) log rfi 



has a minimum point which is stable under C°-perturbations. Therefore, there 
exists a point (d e ,t 6 ) such that V(d,t) Je(se, d e , t e ) = 0. Thus, the function J e has 
a critical point and the claim follows from Proposition 12.31 □ 
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Proof of Theorem \l-4\ In this case 1 = 2 and function ^ defined in (fT8|) reduces to 

*(d,t) = C4<o(fo),l'(&))(*l+*2) 

+ C5< , (&) {(*)"" 2 + (*)"- 2 + 2( ^ 2 )^ 

- c 6 a(^ )(logc?i + logd 2 )- 

It is easy to verify that it has minimum point which is stable under C°- 
perturbations. Therefore, from Proposition ^. 51 we deduce that, if e is small enough, 
the function J e has a critical point. Now the claim follows from Proposition l2.3l □ 

Remark 2.6. The symmetry assumption (<z2) allows to overcome some technical 
difficulties which arise when looking for a solution whose bubbles collapse to the 
same point. Indeed, the problem arises when we study the reduced energy and we 
have to compute the contribution of each peak and the interaction among the peaks. 
The contribution of each peak is clear: it is given by the distance from the peak 
to the boundary as in (|64[) and by the value of the function a at the projection 
of the peak onto the boundary as in On the other hand, to compute the 

interaction among the peaks (see (|65p ) it is important to compare the geodesic 
distance d(si, Sj) between the projections of the peaks onto the boundary with the 
distance \r]iv(si) — T)jv{sj)\ between the normal components of the peaks. To have 
a good expansion the distance d(si,Sj) should be negligible with respect to the 
distance \r]iv(si) — rjjv(sj)\. But then, in order to find a criticality in the points Sj, 
we need to go further in the expansion and computations become too tedious. If 
the domain O and the function a are symmetric, we can overcome this difficulty 
just by assuming that the peaks satisfy (ITT1) . so that d(si, Sj) — 0. In this case the 
interaction among the peaks is clear and it is given in terms of the Green function 
of the Laplace operator on the half-space (see (|5"S1 ), 



Appendix A. Boundary estimates of the Green function 

In this section we establish the technical estimates we used in the previous part. 
We denote by G(x, y) the Green function of the Laplacian with Dirichlet boundary 
condition and by H(x, y) its regular part, i.e. 

G ( x >v) = —, — tt\ — i r^2 ~ H ( x >v)> 

n(n — 2)ui n \x — y\ n 1 

where u> n is the volume of the unit ball in K™. 

First of all, we need an accurate estimate of H(x,y) when the points x and y 
are close to the boundary. Let us introduce some notation. For r\ > we write 
n v '■= {x € SI : dist(a;, 9f2) < rj}. We fix rj small enough so that the orthogonal 
projection p : f^r/ — > 90 onto the boundary is well defined, i.e. so that for each 
x G Q,2r) there is a unique point p(x) £ 90 with dist(x,90) = \p(x) — x\. Set 
d x := dist(x, 90), p x := p{x), and v x := v{x), where as before v{x) denotes the 
inward normal to 90 at x. For x € 02,, we define x := p x — d x v x = x — 2d x v x . 
Thus, x is the reflection of x on 90. 
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Lemma A.l. There exists C > such that 

1 



Cd, 



\x-y 
C 



n-2 



(19) H(x,y)-, — 

(20) K^-^^J -|x-*|- 
for all x € Q v and y G f2. In particular, there exists C > smc/i £/ia£ 

C 



F - y\ 



c 



(21) 0<ff(x,y)< 
and 

(22) |V K ff(x,y)|< 

F-yr 

Proof. For convenience we set 

X(z,2/) := H(x,y) 



x g f2jj, y G f2 



x,y € fi. 



F - y| 



for x G f2„ and y G f2. Note that there is c > 0, only dependent on n and r;, such 
that |x — £| < c\x — £| if x G ri„ and £ G -B(x, d x /2). If moreover y G fi, then 



cd x /2 

If-vl 



< 1 + c, 



(23) ie - 2/1 - 

since y G f2 and dist(£, O) > d x /2. 

The proof of dT9j) is analogous to the proof of Eq. (2.7) in [4], with obvious small 
changes. Similarly, slight modifications of the proof of Eq. (2.8) in .4] yield 

for all x G ri„ and y £ fl. Fix x, y, take r := 2^/n and set 
Q:=UGM" I |a;-£|oo<<Wr}. 
Note that if £ G Q then £ G £?(x, d x /2) and therefore 
(25) (4/2 < di < 3d x /2. 

Hence we obtain for i G {1, 2, . . . , n} 



TTt d 

\d Xi x(x,y)\ < — sup |x(£,y)| + ^sup|A c x(£,y)| 



< C [ sup 



sup ■ 



d x 



< Csup 



iedQ d x \£ - y|»- 2 ?eQ d € |£ - y|«- 2 ^ 
1 



< 



«6Qie-y|"- 2 
c 



\x-y 



n-2 



by [H Eq. (3.15)] 
by HU) and 
by ([251) 
by 121D. 



Summing up this inequality over i gives (|20p. 

To prove (|22[> . note first that there is C > such that 

(26) |V x i2"(x,?/)| < C if x G Q\il v , y G 0. 
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The case x £ £l v relies on the estimate (|20|). Note that there is C > such that 

\% - y\ 



(27) 



F - 2/1 



> C for all a; € 0^, y £ SI. 



This implies that the term on the right of (f2"Uf is estimated by a constant multiple 
of 1/| a; — y\ n ~ 2 if X G and y € O. In view of (|26l) it therefore remains to show 
that 



(28) 



1 



< 



C 



f - yl 



x £ Sl ?; , y G SI 



F - 2/1 

for some constant C > 0. 

Writing ft for d/dxi we calculate as in [3] for any i 6 {l,2,...,n}: 

1 2-n 



(29) 



ft 



^(aij - yj)di. 



3=1 



Since x := x — 1d x v x , we find 
Using this representation in (|29|) yields 



ft 



1 



f - y 



n-2 



< 



c 



\x-y 



—[{I + d x \div x \). 



By our choice of ij we have \d x \ < r\ and \diV x \ < C for all x £ In view of (f2"7| 
we obtain (|28| and finish the proof. □ 



Here and in the remaining appendices we employ the notation 

,1/9 

\u\A,q 

for measurable A C K™ and ge [1, oo]. If A = SI we omit it from the notation. 

Lemma A. 2. Let 6,81,82 £ (0,1] and £,£1,^2 € Sl^. Xe< £ be the reflection point 
of £ respect to dCl. There exists c > swc/i £/iai 

(30) 

and 



< PU 5>t {x) < u s , e (x) 



< U 5ti (x) ~ PUs,(,{x) < a n 8—H{x,() < c- 



(31) ~_~„^~, --«, S v-, „ v »,^^„ |a; _^ n _ 2 
/or all x £ Q. Moreover 

R s ,z(x) ■■= p Us.d x ) - Us,d x ) + a n8^ H{x,t) 

satisfies 

(32) |#5,e|fi,oo = o 
Finally, there is /3 > smc/i £/ia£ 



dist(£,<9Sl)' 



(33) 
(34) 



<5, 
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as 6,61,62 — > 0, independently of £, £1, and £2- 



n 



71-2 



Proof. Estimates ([30]) , (|3lj) , and (|32|) follow easily from the maximum principle and 
Lemma IA.1I 
Note first that 

(35) \U B , c \ q = O (6*-^) if Q > 
and 

(36) Itfjf^O if 9 > 1. 

as (5 — >• 0, independently of £. 
Recall that 



2+2 
i-2 . 



and note that 

(38) V,, < ( 



\x — y\ n 2 



\x-y\ n -i~ 



By ([3"7|). (|2"2l . and ([55]). to show (J33) it; suffices to prove 

71-2 

(39) L L ^•"ff^.w^ - (|) * o 

n+2 

For simplicity, set V := Ug~^ and g(x) := Set M := diam(fi). Pick 

/ n(n — 1) n 
r 6 v (n-l) 2 + f ~1, 

and note that then r > 1 and r' > n, where r' denotes the conjugate exponent of 
r. Since ^7 + i + l = 2it follows as in the proof of [14l Theorem 4.2] that 



u s 2 4 2 (x)g(x ~ y)v(y)dydx < \Us 2 ,£ 3 \ r >\g\B(p,M),r\ v \i 

by (|35p and (|36|) . Here we have used that |<7|s(o, M),r is finite since r < n/(n — 1). 
On the other hand, r > n(n — 1)/((ti — l) 2 + 1) implies that 

/ 1\ 71-2 
71 l-~ = - + P 

\ r ) n — 1 



for some /3 > 0, proving (IMl) and hence 

To prove (|34[) we proceed similarly. This time we pick 

2n \ 2n(n - 1) 



n + 4 ?i 2 + 2?i - 4 



s € max < 1, > 

V I « + 4J 

and define r by 

11 n + 2 

(40) _ + _ = ! + 

r s 2n 
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Some basic calculations reveal that s is well defined and that 

n 



(41) re 1, 

n — 1 / 

Similarly to the proof of [TJ] Theorem 4.2], taking into account the Remark (2) 
following the statement of that theorem, we obtain 

\VPU 5tC \^n_ < \g\ B(0tM) , r \V\ a = 0(5^-^). 

Again we have used that r < n/(n — 1) implies that the r-norm of g in the ball of 
radius M is finite. Since s < 2n(n — l)/(n 2 + 2n — 4), there is (3 > such that 

n n + 2 n — 2 
7 2~ = 2(7i-l) 
proving ([M]). □ 

Appendix B. An estimate of the error 
To simplify notation, from now on we write 

§i '■= <^,ej Ci '■= £i,£j U t \= 

Next, we estimate the error term defined in (Ti~5j) . 
Lemma B.l. It holds true for some a > that 

Ps,d,t|| - O (e^ +a 



Proof. We estimate -R s ,d,t in case (fit)]) . The estimate in case ([7]) is easier and can 
be obtained after minor modifications of this argument. 
From the definition of i* we deduce that 

||i? s ,d,t|| = O (|-div(o(as)VV. jd ,t) - a(x)f t {V a , d , t )\^j 

= O (|-VaW Sld ,t - a(a;)AV; >d ,t - o(ar)/ e (V 8 , d .t)U 



(42) 



a(x) [f{Ui)-f(PUm&, 



O (£\VaVPUi\sA 

+ 0\\a(x)[J2f(PU i )-f(J^PU i )]\_^ ] j 

+ 0(\a(x) [/(T4, dlt )-/ e (V;, d ,t)]|^ 
: h + h + h + h- 



To estimate I\ recall that <5j = O ^e»- 2 J on compact subsets of A. By 
get, for some a > 0, 

(43) \VaVPUi\^_ = O (^ +a ) ■ 
Let us estimate I?,- By (|3T|) for some cr > we obtain 

(44) |a [/(C/,) -/(P^)]|_^ 

= o (\ur 2 (Pu t -u t )\^)+o (\\pui - c^r 1 !^) = o (^ +a ) 
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because by (|3T]) (using also ((48j) with q = (n+ 2)/4) 
(45) \\PU i -U i \*- 1 \_ s ^ = \PU i -U i \ p £ 

n-2 

1 



n + 2 

Si* o 



\x - & 



p-V 



OS 



and by Holder's inequality for some a > (using also (|47p and (14"51) with q 
when n < 6 or q ~ (n + 2)/8 when n > 7) 



IP-2 



(46) 



01 



o 



___2nq__ 

(a-l)(n+2) ■ 



if n > 7 



if n < 6, 



with 



(47) |[/,f- 2 



(„_2)(„ + 2) 



'o{sj) 



o s 



if n > 7 and 1 < q < 2±2, 
if n < 6 and g > 



and 
(48) 



F ^ & 



lra-2 



= o 



2+2\ 

2q 



if n > 6 and q > 1 or n < 5 and 1 < < 



n + 2 
6 — n 



Let us estimate I-i. We set 



(49) 

We have 



16 - 6 



(50) 



a(x) 



\ i ^ i 

+ (E /( p ^)-/(E pu; ' 



n\UiB(ii, v ), 
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because 



(51) 



n\u 4 s(^,ij),^ 



\B(e<,77),^ 



and if j ^ i 
(52) 



Moreover 



+ o (lira-* - c^r 1 ^,,),^) + ofEl^T 1 !^).^) 

and the first term is estimated in (j46)) . the third term is estimated in (|45| . the 
fourth term is estimated in (|5"2"T) . The second term is estimated using (|47|) and (|48|) 
(with q ^ 1 when n < 6 or q ~ (n + 2)/8 when n > 7) as follows 



(53) iC/r'^lste,,),^ 



-°l,««> r W> (n_2)(n + 2) 



o 



1 ['* 



1 ['* 



n— 2— (T^ 



if n > 7 



if n < 6, 



for some cr > 0. 

Arguing exactly as in Proposition 2 of [23], we can estimate the last term I4 by 



(54) 



\a{x) [f (V s , d ,t) - fe (V., dlt )] U = O (e| In e| 

n + 2 



□ 



Appendix C. An estimate of the energy 

It is standard to prove that 

J e (s,d,t) = J e (V s ,d,t) + h.o.t. 

(see for example [3] or [2]), so the problem reduces to estimating the leading term 
Je (^4,d,t) ■ We will estimate the leading term in case (fTOj) . because the expansion 
of the leading term in case ([7]) is easier and can be deduced from that. We also 
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assume I = 2, because with some minor modifications we treat the general case. 
Therefore, the estimate will be a direct consequence of Lemma (|C.3|) and Lemma 

(EH. 

For future reference we define the constants 

(55) 71 = < J ^—L^- dy, 

(56) 72 = < / ; ^^+2 dy, 



( 57 ) 7s = < / „ , , , 2 > n log - ; ^^ dy. 



(1 + M 2 )" 

(i + |y| 2 )" og (t + \ y \*y- 



We start with the following key estimates. 
Lemma C.l. The following estimate holds true: 

(58) J a{x)Uldx^ ll a{s 1 ) + {\/a{s 1 ),v{s l )) ll t 1 e + 0{e 2 ) . 

Here r\ is choosen as in (|49|) . 
Proof. We split the left-hand side as 

(59) J a{x)Ufdx = J a{ Sl )Ufdx+ J (a(x) - afa)) Uf dx . 
b{£i,v) BfaM) Bfa,n) 

We deduce 

(60) J afa)Ufdx = 7l a(£ ) + O (^j . 

By the mean value theorem we get 

a{8\y + £i) - a(£i) = a(5iy + r}iu(si) + si) - a(£o) 

= (Va(si), i/( Sl )> 77! + <5i (Va(si), y) + R(y), 
where R satisfies the uniform estimate 

(62) \R{y)\<c(8l\y\ 2 + 8^^+111) for any y G B(0, rj/8i). 

Therefore we conclude 



(63) J (a(x) -a(^))Ufdx 

1 

(i+Ty 



a(£y + T7ii/(si) + s x ) - a(si)] jr-^dj 



cr 

B(0,»j/<5i) 

a£ y K Vo(ai), v( Sl )) m + ft < Vo(ai), y) + R(y)} + j 2 dy 

3(0,77/5!) 
(Va( Sl ),Ksi))7im+0(r; 2 ). 

□ 
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Lemma C.2. The following estimates hold true: 
(64) 



J a{x)U{- 1 (PUt - U t ) dx = - l2 a( Sl )e (J^j + O (e 1+ °) 



and 
(65) 



a{x)Ul~ l PU 2 dx 

' O (e 1+ °) 
= < l2a{io)e{did 2 )~ x 



|ti-i 2 |"- 2 |ii+t 2 |"- 2 



ifsi ^ s 2 , 



0(e l+a ) if Sl = S2 =i Q , 



for some a > 0. Here r\ is choosen as in (|49[) . 

Proof. First we prove (|M|) . By Lemma \A. II and Lemma [A. 2 1 we get 



(66) / a{x)Ur l {PUi - U{] dx 



a{x)Ul- x {-a n S^ 2 H(x,b) + 



dx 



3(0,77/5! 



a{hy + £i)H(6iy + 6,6)- , — , 1+2 dy 

(1 + \y\ 2 ) 2 



B(0,?j/«i 



a(5iy + 6) 



1 



1 



|^iy + £i-£i| n - 2 (l + M 2 )^ 

n—2 \ / / r \ n- 

<5l \ 



-dy 



»7i 



v 1 f . 



2771 



a(si) 



O 



1 



(i + M 2 )- 

n-l\ 



O 



n-2 



»7i 



because 

l^iy + 6 - 61 = |*i« + 2»jii/(*i)| > 2 m - |<5x2/| > »?i for any y £ B^/Sx). 
and by mean value theorem 



a(£iy + £i) = o(si) + O (?7i) and 



1 



■O 



i^+a-6i n - 2 (2??i)^ 2 v<-\ 
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Next, we prove ([65]). By Lemma FO 
(67) J a(x)Uf~ 1 PU 2 dx 

a(a;)t/f _1 (u 2 - a n 5 2 2 H(x,£ 2 ) +-R<5 2 ,£ 2 ) dx 



B(ei,f)) 



(i + M 2 )^ 

B(<W«l) 

" + -61")-* + 



(i + M 2 )^ 

B(0,77/5i) 

1 1 



i-2 



+ (^) + o(^,) + o(<^f). 

because for any y G B(0,r]/8i) we have 

|<5i2/ + 6 - &| = l<5i£/ + (»7i + ?72M£o)l >Vi + V2- \5ty\ > r) 

\6iy + Si-b\>\ti-S2\-\Siv\>v 
and by mean value theorem a(d±y + £1) = a(£o) + (?7i) and 

1 1 



lr/i - m\ n 2 |r?i + m\ n 2 V »7 



n-2 



i i ,, fSi\y\ + s- 2 



,Tl-l 



□ 
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Lemma C.3. The following estimate holds true: 

Jo (14,d,t) = ^ [27ia(Co) + 7i (Va(£ ), e(ti + t 2 )] + ^7 2 a(£o)x 

di\"~ 2 ( d 2 \ n - 2 l l 

+ 0(e 1+ff ) 

/or some cr > 0. 
Proof. 

(68) JoG4, d ,t) = \J a{x)\VV s . A . t \ 2 dx~- J a{x)\V s ^ t \ p dx 

n n 

We estimate the first term at the R.H.S. of (1681). We write 



(69) / a(x)|W d , t | 2 cfe 

= J a(x)\VPUi\ 2 dx + J a(x)\VPU 2 \ 2 dx-2 J a(x)V PU X V PU 2 dx 

n n n 

Let us estimate the first term in (|69[) . The estimate of the second term is similar. 
Let us choose 77 as in (|49"j). We get 

y a{x)\\7PUx\ 2 dx = - J div {a{x)\7 PU X ) PU x dx 
n n 

= - J (a(x)APUx) PUxdx - / (Va, VPC/i) PJ/icfo 

SI 12 

( ?0 ) = y a{x)Ul~ l PUxdx- J (Va,VPUi) PUxdx 



a(x)Uf~ x PUxdx + I a^ur 1 PUxdx 

-J (Va,VP(/i) PUxdx 
n 

By we deduce for some /3, cr > 

(71) y (Va, VPf/i) PUxdx <C J \VPUx\PUxdx = O (s}^ +fi J = O (e 1+CT ) . 
By Lemma IX2I we also deduce 

(72) y a(x)Uf~ 1 PUi<kc = ((j- 

n\B(.t lt ri) 



BOUNDARY CLUSTERED LAYERS NEAR THE HIGHER CRITICAL EXPONENTS 21 



and 
(73) 



J aixjUr 1 PUidx = J a{x)Uldx+ J a{x)Uf~ l {PUi - U x ) dx. 



The first term is estimated in Lemma IC.ll and the second term is estimated in 
of Lemma IC.2I 

It remains only to estimate the last term in (|69p. 

J a{x)V 'PU X V , PU 2 dx = - J div (aVP£/i) PU 2 dx 

( 74 ) = - J {aAPU 1 )PU 2 dx - J (Va, VPU X ) PU 2 dx 

a si 

= j a(x)Uf' 1 PU 2 dx- J {Va,VPUi)PU 2 dx. 



We have 
(75) 

and 



J a(a;)t/f~ 1 PU 2 dx = J h / 



(76) / a(x)f7f _1 PC/ 2 dx 
n\B(ei,r,) 







ti + 2 n-2 


/ 




B(«i,»?) 


n+2 n-2 




77™ 


/ 

\B(0,1) 











1 



-dx 



\y\ n+2 \y + ^K 2 



dy 



The first term in (|75|) is estimated in (|65|) of Lemma IC.21 
Finally, as in the proof of ([7T|) , from ([33j) we obtain 

(77) / (Va, VPE/ 2 ) PJJida; = O (e 1+CT ) , 

since < C\ < 8 2 /8i < C 2 on compact subsets of A. 
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We estimate the second term at the R.H.S. of (l68l). We write 



f a{x)\V d , t \ p dx = J a{x)\PU x - PU 2 \ p dx 
n n 

(78) = / a(x) (\PUi - PU 2 \ P - \Ui\* - \U 2 \ P ) dx 

a(x)(\U 1 \P + \U 2 \P)dx. 



The last two terms in (f78|) are estimated in Lemma [C. II Let us choose r\ as in 
49}. 

We split the first integral as 



(79) / a(x) QPUi - PU 2 \ P - \Ux\ p - \U 2 \ P ) dx 

... + / ... + / 

B(tuv) B(&,ri) 0\(B(Ci : ,,)Ui3(C2,.))) 

From Lemma lA.2l we deduce 

(80) j a(x)(\PUx-PU 2 \ p -\Ux\ p -\U 2 \ p )dx 



I 

= O 



\p\(BKl,tj)UB(6,»7)) 



We now estimate the integral over B(^i,r]). 
(81) J a(x)(\PUx-PU 2 \ p -\Ux\ p -\U 2 \ p )dx 

= p j aixjUf' 1 (PUx -U x - PU 2 ) dx 



2 

B(Ci,v) 

a(x)\U 2 \ p dx 

B(ii,r,) 



J a(x)\Ux + 9 (PUx -Ux- PU 2 ) \ p ~ 2 (PUx -Ux- PU 2 f dx 



P 



J a(x)J7f _1 (PUx - Ux - PU 2 ) dx + I, 
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where / is defined and estimated as 
(82) 

[p - l)p 



I := 



J a(x)\U! + 9 (PUi - U± - PU 2 ) | p ~ 2 {PUi -U x - PU 2 f dx 



a(x)\U 2 \ p dx 



= o 



O 



J Ur 2 (PUi - Uif dx\+o\ J Uf- 2 Uldx 
J \ PUl -uA + o( J \U 2 \Pdx 

\B(£i,v) J \B(ii,v) 

P-2 



= O [\ur z (PUi ~ U X )\ B{ ^^\PU X - u x \ B{i ^^ 

+o(ic/r 2 f/ 2 i B(Cl ,,),^if/ 2 i s(6 ,, )) ^) 



for some a > 0, because of estimates (|35|) . (05]), (1521 and (155)) . 

The first term in (f5Tj) is estimated in (|5I|) and (1551) of Lemma IC.2 
Finally, we estimate the integral over B(^ 27 -q) 
(83) r 

a{x)(\PU 1 -PU 2 \P-\U 1 \P-\U 2 \P)dx 



-p / a(s)E/f _1 (-PU 2 + U 2 + PU X ) dx 



B(ti2,v) 

/ a(^)|{7 1 +0(-PC/ 2 + t/ 2 +P;7 1 )r 2 (- J P;7 2 + C/ 2 + PC/ 1 ) 2 d a; 

B(?2,'f) 

a(x)|t/i| p da; 
B(6.n) 

= p / a{x)U^ [PU 2 - C/ 2 - PUx) dx + J, 



where J is estimated exactly as in (|82p . while the first term in ([83[) is estimated in 
(|63j) and (|65j of Lemma 



We collect all the previous estimates and we get the claim. 



□ 
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Lemma C.4. The following estimate holds true: 

(84) — [ a(x)\V s , d , t \^dx = - f a(x)\V s , d , t \vdx 
p-e J P J 



4 J a(x)\V s , d ,t\ p dx-- J aix^V^-HoglVs^tldx 



o(e) 



r / s , / Ni I 71 71 «n 73 , 

= [a(si) + a{s 2 )\ I — I e 



p p 
n — 2 

+ -7: — 7i Ksi) log Si + a(s 2 ) log S 2 ] e + o(e). 
Proof. We argue exactly as in the proof of Lemma 3.2 of [5]. 



□ 
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